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Abstract. HELAC-PHEGAS is a FORTRAN based package that is able to compute automatically and efficiently tree- 
order helicity amplitudes and cross sections for arbitrary scattering processes within the standard electroweak theory and 
QCD. The algorithm for the amplitude computation, HELAC, exploits the virtues of the Dyson-Schwinger equations. 
The phase-space generation algorithm, PHEGAS, constructs all possible kinematical mappings dictated by the amplitude 
under consideration. Combined with mutichannel self-optimized Monte Carlo integration it results to efficient cross 
section evaluation. 



INTRODUCTION 

The need for efficient algorithms to calculate helicity 
amplitudes and cross sections for any process, in an au- 
tomatic way, has been well recognized long time ago. 
Up to now, algorithms that efficiently combine helic- 
ity amplitude computation and phase-space integration 
have been prooven successful for specific processes, like 
for instance four-fermion (1) production in e+e^ col- 
lisions. On the other hand general-purpose computa- 
tional packages Uke CompHEP (2) and GRACE (3) do not 
provide automatic efficient phase-space integration algo- 
rithms. Moreover the vast majority of the automatized 
helicity amplitude computationial algorithms, like for in- 
stance MadGraph (4), have been based on the Feynman 
graph representation of the amplitude which severely re- 
stricts their ability to deal with multiparticle scattering 
processes. 

In this article we report on some developments that 
have lead to the construction of two programs, HELAC (5) 
and PHEGAS (6), that allow for an efficient and automatic 
evaluation of cross sections for arbitrary scattering pro- 
cesses. 



therefore as n I where n is the number of particles involved 
in the scattering process. 

An alternative^ to the Feynman graph representation 
is provided by the Dyson-Schwinger approach. Dyson- 
Schwinger equations express recursively the n-point 
Green's functions in terms of the 1 — ,2— ,...,(« — 1)- 
point functions. For instance in QED these equations can 
be written as follows: 





+ ^ {ig)n>iii,iP2)jMPiHPuP2) 

P=P[+P2 



where 
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HELAC 

The traditional representation of the scattering ampli- 
tude in terms of Feynman graphs results to a computa- 
tional cost that grows like the number of those graphs. 



describes a generic n-point Green's function with respec- 
tively one outgoing photon, fermion or antifermion leg 
carrying momentum P. n^,v stands for the boson prop- 
agator and e takes into account the sign due to fermion 
antisymmetrization. 

In order to actually solve these recursive equations it is 
convenient to use a binary representation of the momenta 
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See also references (7) and (8). 



involved (8). For a process involving n external particles, 
all momenta appearing in the computation, Pi^, 

iei 

where / C { 1 , . . . , n}, can be assigned a binary vector in = 
(mi , . . . ,m„), where its components take the values or 
1, in such a way that 

(=1 

Moreover this binary vector can be uniquely represented 
by the integer 

n 

m = ^I'^^nii 

i=l 

and therefore all subamplitudes can be labeled accord- 
ingly,i.e. 

A very convenient ordering of integers in binary repre- 
sentation relies on the notion of level /, defined simply 
as 

n 

I = Y^m . 

1=1 

As it is easily seen all external momenta are of level 1, 
whereas the total amplitude corresponds to the unique 
level n integer 2""^ This ordering dictates the natu- 
ral path of the computation; starting with level- 1 sub- 
amplitudes, we compute the level-2 ones using the 
Dyson-Schwinger equations and so on up to the level n 
one which is the full ampUtude. For the spinor wave func- 
tions as well as for the Dirac matrices, we have chosen the 
4-dimenional chiral representation which results to parti- 
curarly simple expressions. All electroweak vertices in 
both the Feyrmian and the Unitary gauge have been in- 
cluded. 

The computational cost of HELAC grows like ~ 3", 
which essentially counts the steps used to solve the recur- 
sive equations. Obviously for large n there is a tremen- 
dous saving of computational time, compared to the n! 
growth of the Feynman graph approach. 

For QCD amplitudes colour representation and sum- 
mation plays an important role. Let 1 . . . n denote the 
colour labels of quarks and o,(l) . . .o,(«) denote the 
colour labels of antiquarks, with = \ ...n\ being 

a permutation of {1 . . .n}. The colour factor is given ob- 
viously by 

Cl = 5la,(l)^2ai(2) • • -SnaKn) 
Moreover the colour matrix, defined as 

colours 



with the summation running over all colours, \ ...Nc, has 
a very simple representation 

where m{Gi,Cj) — 1 counts how many elements of the 
permutations o, and Oj are common. In order to extent 
this colour representation to QCD amphtudes we have 
just to consider gluons as being quark-antiquark pairs 
and assign to them two colour labels (/,0,). The colour 
factor and the colour matrix still has exactly the same 
form. The only thing one has to consider is to rewrite 
the known Feynman rules of QCD in a slight different 
way. It is worthwhile to note that exact colour summa- 
tion is efficient as far as the number of equivalent gluons 
is smaller than O (5 — 6). For multicolour processes other 
approaches have to be considered (9, 10). 

The programme is also incorporating the possibiUty 
to use an extended precision by exploiting the virtues of 
FORTRAN90. The user can easily switch to a quadruple 
precision or to an even higher, user-defined precision by 
using the multi-precision library (11) included in HELAC. 
In this way, a straightforward computation of cross sec- 
tions for processes Uke e~e~^ — > e~e~^e~e'^ without any 
cut is reUably performed (12). 



P HE GAS 

The study of multi-particle processes, like for in- 
stance four-fermion production in e+e^, requires efficient 
phase-space Monte Carlo generators. The reason is that 
the squared amplitude, being a comphcated function of 
the kinemtaical variables, exhibits strong variations in 
specific regions and/or directions of the phase space, low- 
ering in a substantial way the speed and the efficiency of 
the Monte Carlo integration. A well known way out of 
this problem rehes on algorithms characterized by two 
main ingredients: 

1. The construction of appropriate mappings of the 
phase space parametrization in such a way that the 
main variation of the integrand can be described by 
a set of almost uncorrelated variables, and 

2. A self-adaptation procedure that reshapes the gener- 
ated phase-space density in order to be as much as 
possible close to the integrand. 

In order to construct appropriate mappings we note that 
the integrand, i.e. the squared amplitude, has a well- 
defined representation in terms of Feynman diagrams. It 
is therefore natural to associate to each Feynman diagram 
a phase-space mapping that parametrizes the leading vari- 
ation coming from it. To be more specific the contribu- 
tion of tree-order Feynman diagrams to the full amplitude 



Table 1. Results for several processes using HELAC-PHEGAS. In the second column the number of Feynman 
graphs and in parenthesis the number of steps required to solve the recursiveDyson-Schwinger equations are 
given. 



Final states 


Number of FG(DS) 


CGeV) 


Cross section ffb) 


e~ ud s cy 


90(74) 


200 


199.75 (16) 


e-e+ ^e- VeH+V^y 


108(100) 


200 


29.309 (25) 


e~e+ -^iT Vf,u dyj 


587(210) 


500 


1.730 (58) 


e~e~^ — > /i~ Vju M dc c 


209(102) 


500 


0.1783 (20) 


e~e+ — » /i~ M d c c y 


2142(339) 


500 


0.02451 (65) 


gg-*bbbbW-W+ 


960(380) 


500 


4.716(24) 



can be factorized in terms of propagators, vertex factors 
and external wave functions. In general, the main source 
of variation comes from the propagator factors and there- 
fore our aim is to construct a mapping that expresses the 
phase-space density in terms of the kinematical invariants 
that appear in these propagator factors. Since in prin- 
ciple we need as many mappings as Feynman diagrams 
for the process under consideration, we have to appropri- 
ately combine them in order to produce the global phase- 
space density. A simple and well studied solution to this 
problem was suggested some time ago in reference (13). 
It should be mentioned however that other self-adapting 
approaches can be used as well (14). It is important to 
note that although by using Feynman graphs to construct 
phase-space mappings we face the original n! computa- 
tional cost growth problem, the self-optimization cures to 
a certain extent this by selecting only the few mappings 
that dominate the phase-space density. For alternative ap- 
proaches we refer to (15). 

In order to describe the construction of the phase- 
space mappings, let us consider a typical process in which 
two incoming particles produce n outgoing ones. The 
phase space, d^n{P = qi+q2',Pi--- ,Pn), can be decom- 
posed as follows 

d^n=\fl^]d^^{P;Qi,...,Q„,) 



,i=i 



where the subsets {n,?"2, • • • ,''ni} up to {si,i2, • • • 
represent an arbitrary partition of {p\,pi, ■ ■ ■ ,Pn}- The 
above equation can be generalized recursively resulting 
in an arbitrary decomposition of d^„. Feynman graphs 
can be seen as a reaUzation of such a decomposition, this 
latter being identified with a sequence of vertices of the 
graph. There two possible cases for 2 ^ n scattering. 



First, all outgoing momenta involved in the vertex are 
time- like. 




d^„ = ...^^dMQ^QuQi)... 

dQ\ dQl 'k'/'-iQ^QlQl) 
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with A,(x, >', z) = + — 2xy — 2xz — 2yz, and second 

when one of them is space-hke. 



d<^n 



dQl dQl 

2% 2% 
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with 



= {Qi-qif 

= ml + Ql-^{Q' + Ql-Ql) 



Q 



Q 



^2 1 cos 9 



and (£2,^2) being the incoming momentum q2 in the 
rest frame of Q. The appropriate sequence of vertices. 



{Vi,V2, ■ ■ • , Vjt} can be chosen in such a way that a re- 
cursive construction of the phase space is reaUzed. For 
instance Vi should contain at least one incoming particle 
whose momentum is known. The rest of the sequence is 
chosen recursively: vertex Vj is characterized by an in- 
coming momentum Q which has already been generated 
in one of the {Vi, . . . 

Following the above described algorithm we end up 
with an expression for the phase-space density, 

where s,- and tj refer to the kinematical invariants enter- 
ing the propagator factors of the graph and and cos 6; 
represent center-of-mass angles needed to complete the 
phase space parametrization. It is now straightforward 
to generate Si and tj with probability densities (Pi{si) and 
^j{tj) that are automatically chosen accordingly to the 
nature of the propagating particle. 

Results, demonstrating the ability of 
P HE GAS -HE LAC to deal with multiparticle processes, 
are presented in table 1 (6, 16). 



SUMMARY AND OUTLOOK 

PHEGAS-HELAC offers a framework for high-energy 
phenomenology. It provides all necessary and suffi- 
cient tools for efficient, reliable and automatic compu- 
tation of helicity amplitudes and cross section. The Stan- 
dard Model, including QCD, has been fully incorporated. 
Higher-order corrections are in principle tractable within 
the framework of Dyson-Schwinger equations and work 
is in progress in order to include electroweak corrections 
as described in reference (17). New physics interactions 
and models, icluding the Minimal Supersymmetric Stan- 
dard Model and the trilinear gauge couplins will be con- 
sidered in the near future. 
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